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Abstract 

We study three-body charmless baryonic B decays of H ^ AAP with P = n and K in the 
standard model. We find that the branching ratios of the K modes are about one order of magnitude 
larger than those of the corresponding vr modes unlike the cases of B ^ ppP. Explicitly, we obtain 
that Br{B- AAR-) = (2.8 ± 0.2) x 10"® and Br{B^ AAi?®) = (2.5 ± 0.3) x 10"®. The 
former agrees well with the BELLE experimental measurement of (2.91^'^Q'yQ ± 0.38) x 10“®, while 
the latter should be seen at the ongoing B factories soon. 

PACS numbers: 
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There have been lots of attentions recently on charmless three-body baryonic B dec^s 
due to the several new experimental measurements by BELLE and BABAR Q, Q, S, 0, . 

It has been realized that baryonic B decays could happen since the R-meson mass can be 
heavier than the invariant baryon masses. In particular, the reduced energy release can make 




the decays as signihcant as the two-body mesonic decay modes m . As the main characteristic 
of the three-body baryonic decays, the baryon pair threshold effect j6|, l7| results in the decay 
modes being accessible at the current B factories. Indeed, some of the modes have been seen 
and limited recently by the BELLE and BABAR collaborations, with the data given by 
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(3.06;[|:^^ ± 0.37) X 10“® (BELLE) 

± 0.14) X 10"® (BELLE) 
(5.3l^:3^ ± 0.58) X 10“® (BELLE) 

(6.7 ± 0.9 ± 0.6) X 10"® (BABAR) j3| , 
(3.27l0;®2 ± 0.39) x 10“® (BELLE) 
(2.16l[(|^±0.20) X 10"® (BELLE) 
(2.9ll[J;?J( ± 0.38) x 10"® (BELLE) 


(1) 


and 


IP. 


Br{B- AAtt-) < 2.8 x 10"® (BELLE), (2) 

respectively. Note that the three-body decays in Eq. o are much larger than the cor¬ 
responding two-body baryonic modes, for which only upper bounds of 0(10“^) have been 
reported js|. 

There are mainly two kinds of approaches to study the baryonic B decays in the literature. 
One is the pole model, presented in Refs. 9|, [lOl, [llj, where the intermediate particles 
couple dominantly to the final states. The other, proposed in Refs. R llE 13) is based 
on the QCD counting rules 3 3) which deal with the baryonic form factors by power 


expansions. We note t 
and Br{B~ —>■ ppK~) 


Q 


mt global fits in Ref. |3 for Br{B ppn ), Br{B^ ppK^) 


with the QCD counting rules have been performed and consistent 


results with data have been derived. Furthermore, various radiative three-body baryonic B 


decays 


Q 


have been studied. 

In this report, we will concentrate on the three-body charmless baryonic decays of 

B AAP (F = tt and K) (3) 
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FIG. 1: Diagrams for B —> BBP with B = p and A, P = 'ir{K) for q = d{s) and q' = u{d) for 
charged (neutral) modes, and where (a) and (b) represent C and T terms, respectively. 

based on the QCD counting rules. We note that so far there has been no theoretical study 
on the modes in Eq. m in the literature. In particular, we do not know the reason why the 
decay branching ratio of B~ —AAvr” is smaller than that of B~ —> AAK~, whereas the 
corresponding pp modes are comparable, in terms of the data shown in Eqs. and ((21). In 


addition, these c 
violating effects 


ecays are of great interest since they provid e op portunities for probing T 
1^ due to the measurable polarization of A |lf| . 

There are two types of diagrams which contribute to the decays in Eq. dSD under the 
factorization approximation |2|, with Figs. Et and Eb representing typical diagrams of a 
current-produced baryon pair with B P and a current-produced P with B BB, 
named as C and T terms, respectively. The amplitude oi B BBP with B = p or A and 
P = TT oi K in the factorization approximation is given by [3] 


A{B BBP) = C{B BBP) + T{B BBP) , 


( 4 ) 


where [19|, El 

C(P-^ BBP) = ^{Vuiy:M^^\{nu)v-A\Q){P\mv-A\B) 


Gf 


V2 


-VtiM:, 


a 3 (BB|(MM + dd + ss)y_^|0) -h a 4 {BB\{qq)v-A\ 0 ) 


-l-a5(BB|(MM + dd + ss)y+A|0) -l- -a9(BB|(e„MM -f eadd + esss)y_^|0) 


(P|(g6)y_A|P) + E,,E42a6(BB|(gg)5+p|0)(P|(gfe)5-p|P) , 


( 5 ) 
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with q = d{s) for P = 7i{K), while T{B ^ BBP) is decomposed with the P-meson decay 
constant fp induced from the P creation and B BB through scalar and pseudoscalar 
currents as 


T(P ^ BBP) = 


ap{BB\q'b\B) + /3p{BB\q'^5b\B) 


( 6 ) 


with g' = u{d) for charged (neutral) modes, and 


Oip{Pp) — 


VubV*ap — VtbVf* ( a'p ± a' 


2m?p 




(7) 


where — ^k~ — — ^ 2 ? — 0? ^ 7 t~ — ^k~ — d(X 9 / 2 , 

= Og and a"o = —Here, = /jt/a/^, Vij denote the CKM mixing matrix 


a" 


7 t -, k-,ko 

elements j2l|, Oj (i = 1, • ■ ■ , 10) are composed with (effective) Wilson coefficients c. 
1, 2, • • ■ , 10), defined in Refs 
Oj and related by 


>//) 


• QQ, 


= 


and Nc is the effective color number. The coefficients 


„e// e// 

eff I W+1 /• eff I ^i—1 /■ \ 

tti = c/-' + (*=odd) , tti = c/-' + —— {i=even) 


( 8 ) 


To calculate the decay rate, we need to know the hadronic transition matrix elements in 
Eqs. 0 and Q. The B ^ P transition matrix element can be parameterized as 


(P|g7^(l-75)^|H) = 


(Pb + PpT 


m 


B 


nil 


{Pb - Pp)' 


:{Pb-PpY 




m 


B 


ml 


iPB-ppYPY^Yt), 


{Pb - PpY 

where t = (pe + PbY are defined by j2 

Pf^^(O) 


(9) 


pr^t) = 


(7 JP')Y M, 


P'B^P 


it) = 


TYB^P 


_|_ £1^' 


( 0 ) 


(JQlt _|_ (702 ’ 

Ilf 


( 10 ) 


with the input parameter values of Ff^^{0) = P(f^’^(0) = 0.29, an = 0.48, an = 0, 
(Joi = 0.76, ao 2 = 0.28 and My = 5.32 GeV for P —»■ vr and FY^^{0) = P(f^'^(0) = 0.36, 
an = 0.43, a^ = 0, doi = 0.70, cro 2 = 0.27 and My = 5.42 GeV for B —>• K, respectively. 
For those of the baryon pair involving the vector, axial-vector, scalar and pseudoscalar 
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currents in Eq. we have 


(BBIV'^IO) = + +Pb)m|u(pb) , 

me + rng j 

= h(pB)([^i(t) + ^2(i)]7M + — {Pb - Pb)Av{pb) , 

[_ fn-B + uib J 

(BB|7l^|0) = m(pb)(^a(^)7m + {Pb + Pb)Ai5v{Pb) , 

I mB + rriB j 

(BB|^|0) = fs{t)u{pB)v{pB) , 

(BB|P|0) = gp{t)u{pBh5v{PB), (11) 


where = qa^^Qj, S = q^qj and P = qa^qj with qi = u,d and s. We note 

that F 2 {t) alone can not be determined by the present experimental data. However, F 2 {t) 
can be ignored since it acquires one more 1/t than Fi{t) according to the power expansion 


in a perturbative QCD re-analysis j23, 
quark mass limits, we have 


24| . By using equation of motion and adopting zero 


hA{t) = -gA{t) 


{tub + m-Bf 


( 12 ) 


In terms of the SU{3) flavor symmetry, the form factors [Fi(t) + F 2 {t)], gA(t), fs(t) and 
gp{t) in Eq. can be related to another set of form factors Dxid')^ Fx(t') and Sx(i') 

where X = V,A, S and P denote vector, axial-vector, scalar, pseudo-scalar currents, defined 
in Table m respectively. It is noted that the zero value of (pp|(ss)x|0) in the second column 
is due to the OZI suppression rule. As an illustration, we take (AA|(-un + dd + ss)y,A|0) in 


TABLE I: Relations of form factors between different sets of parameterizations, where the form 
factors [Pi(t) -f F 2 {t)], gA{t), fsit) and gp{t) correspond to the notations V, A, S and P in the 
first column, respectively. 


Form Factor (X=V,A,S,P) 

BB = pp 

BB = AA 

(BB|(ttu)x|0) 

Dx{t) + Fx{t) + Sx{t) 

l[Dx{t) -h35x(t)] 

(BB|(dd)x|0) 

Sx{t) 

l[Dx{t) +3Sxit)] 

{BB\{ss)x\0) 

Dx{t) - Fx{t) + Sx{t) = 0 

li-^Dxit) + 3Sx(t)] 

(BB (■uti + dd + ss)x 0) 

2Px(t) + 3Sxit) 

2Px(t) + 3Sx{t) 

{BB\{euUU + eadd + esss)x 0) 

Fx{i) + \Dx{t) 

-iDxit) 
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Eq. (0 and we have 


FiHt) + F 2 ^Ht) = ‘^Dv{t) + 3Svit) , 
9A^it) = 2 E>A(i) + 3S'A(t) , 


( 13 ) 


For Dxit), Fx(t), and Sxit), since they are related to the nucleon magnetic (Sachs) form 
factors we adopt the results for in Ref. j^ . which are extracted from 

experiments. 

The functions of Dxit), Fx{t) and Sx{t) are parameterized as 


Gl,(t) = J2 


2=1 


Xj 


ln( 


^'■0 


'' Sv(i) = E tSt [1 

i=l 0 


i=l 0 


2=1 

in(^) 

L ^'■0 . 

-7 ^ f 

i=l ^ 

In(^) 

L ^'■0 J 

2=1 

In(^) 

L ^'■0 J 

2=1 

L ^'-0 J 

-7 ^ f 

.Fp(t) = E-^ 
2=1 

L ^'-0 J 

—'y 2 _ 

2=1 

M^) 

L ^'■0 J 


-7 


-7 


-7 


(14) 


with 


Dvit) = , Fv{t) = Gl,(t) + iGi(«), Ds(t) = -^n,Gl,(t ), 

Fsit) = n,lG^(() + 1 g;,(()| , Ssit) = n,Sv(t). 


(15) 


The input numbers can be found in Refs. [7|, [12, ll^. Explicitly, we take 7 = 2.148, 
xi = 420.96 GeV^ X 2 = -10485.50 GeV®, X 3 = 106390.97 GeV®, = -433916.61 GeV^°, 
X 5 = 613780.15 GeV^^ yi = 292.62 GeV^ 2/2 = -579.51 GeV®, si = xi - 2 yi, S 2 = 
500 GeV®, di = xi - ||/i, /i = |xi + ||/i, Si = -Xi/3 + 2 yi, = f 2 = -478 GeV®, 
S 2 = 0, di = Hq^yi, fi = Uqixi - ||/i). Si = Ugixi - 2yi), sg = n^Ss, ^2 = /2 = -952 GeV® 
and Aq = 0.3 GeV. We remark that the parameter Uq in Eq. fTH|l corresponds to the 
(ruB — xn-Qi)/{mq — rriq/) term in connecting scalar form factors to vector ones in the case 
of B ^ B'. In R ^ BBP decays, this parameter is not well-defined. However, by taking 
rriqi —> rriq and mg/ —> me, it has been shown that Uq is around 1.3 — 1.4. Here we fix 


Uq ~ 1.4. 
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For the B —> BB transition in Eq. ©, we first discuss the case with B = p. For 

□ 

B —>• pp, one has |7( 


{pp\ub\B ) = iu{pp)[FA + Fp75\v{Pp) , 
(pp|m75 &|5”) = iu{pp)[Fv Fs]v{pp), 


(16) 


where p = Pb — {pp + Pp)- Note that Fs = Fp as shown in Ref. For B^ —> pp, one gets 

{pp\db\B^) = iu{pp)[F^^ 5 + F^P-f5]v{pp) , 

{pp\dj5b\B^) = iu{pp)[FPP ^+Fs^]v{Pp)- (17) 

The form factors in Eqs. (1161) a nd (HZD can be related |3, by the SU(3) symmetry and 
the helicity conservation \14. Il5| and one obtains that 


F^=^illFA + 9Fy), Ff = 1(9F^ + 11F^), =—Fp. 


10 


10 


(18) 


Moreover, the three form factors Fa, Fy and Fp in Eq. (HED can be simply presented by 


Q 


Fap = 




V 




(19) 


fi ' ' 

where Cj (i = A, V, P) are new mrameterized form factors, which are taken to be real. By 
following the approach of Refs. P, the form factors in B~’^ —> AA transitions are given 
by 


p-iAA 

ipAA _ '-^A 


pyAA _ pyAA 

zpAA _ '-'V zpAA _ '-^P 

^P(S) - 


( 20 ) 


where 


cy = i(9C^ + 6CA , cy = i(6C^ + 9Cv) , 0^ = 0. 


( 21 ) 


It is interesting to note that = 0 from Eqs. (j20jl and (j2H)- The reason^ for this result is 
that 7^p(^) correspond to a helicity (chirality) flipped terms in the B —> AA matrix elements 
in the large momentum transfer. It is well known that the spin of A is carried by the s-quark 
component. Hence, such a term requires a chirality flip in the s-quark, which is absent in 
the decay amplitude oi B —>■ AA. 


^ We thank the argument given by the referee. 
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TABLE II: Fits of {Ca^Cv) in units of GeV^. 


Input 

experimental data 

Fit result 

best fit (with la error ) 

Br{B~ pp7T~) 

Br(B0 ^ ppK^) [2] 

Br{B~ — >ppiL“)[2j 

Br[B~ —> Apy) 

Eq. ([U) 

Ca 

Cv 

X^/ndf 

-71.0 ±5.3 

42.0 ±9.2 

1.9 


To obtain these unknown form factors, we use the fif with the experimental data in 


Eq. (^. Here we have neglected the Cp term since it has one more 1/t over Ca and 
Cy as shown in Eq. m- Therefore, we keep the numbers of the degree of freedom (ndf) 
to be 2. In our £t, we also include the uncertainties from the Wolfenstein parameters 2g 


in the CKM matrix. Explicitly, we use A = 0.2200 ± 0.0026, A = 

Te Refs. 

23. 


p = 0.20 ± 0.09 and p = 0.33 ± 0.05 j2 


di 


/A2 = 0.853 ± 0.037, 


20l. l28j | to deal with a* 


We follow 

[i = 1,..., 10) and we take Wilson coefficients from Ref. 

The experimental inputs and the htted results are shown in Table El With the htted 
values in Table m the theoretical predictions with la error for B AAP are given by 


Br{B- AAR-) = (2.8 ± 0.2) x 10"® , 

Rr(R“ ^ AAtt”) = (1.7T0.7) X 10"^ (22) 

which agree very well with the recent BELLE data P| in Eq. ©• We note that the 
spectrum of B~ AAK~ shown in Fig. [2b is consistent with that in Ref. by BELLE. 
Our predicted values in Eq. (|22I) show that Br{B —^ AAvr ) <C Br{B AAK ) unlike 
Br{B~ —> pp7r~) ~ Br{B~ —>■ ppK~) as indicated in Eq. (^. This result can be easily 
understood from the theoretical point of view. As B~ —»■ pp'K~ can not obtain a large 
contribution from C in Eq. © since 02 is color suppressed, its main contribution is from oi 
term in T as seen in Eqs. and whereas B —> ppK is due to the penguin part which 
gives contributions to the terms in both C and T. However, it is not the case for B~ —> AAK~ 
since (AA|(ss)x| 0) in Eq. ® escapes from the OZI suppression. Its contribution is mainly 
from C which is enhanced by with the chiral enhancement, whereas it is suppressed in 
T. For B~ —AAtt”, on the other hand, because of the 02 color suppression in C and 
the small contribution in T, its branching ratio is small. Moreover, to explicitly see the 















FIG. 2: dBr/dm-Q-Q as a function of m^g for (a) B ppP and (b) B —> AAP, where the solid 
and dashed lines for P = K~ and P = 7r~, respectively. 

T-suppression in the AA case, we show the spectra of B~ —> pp'K~ and B~ —> AAtt” in Fig. 
121 We note that the two modes contain the same T-amplitudes in Eq. ©. As shown in 
Fig. El the main contributions to the rates are due to the threshold effect, resulting from 
the form factors of FA,y = CAy/t^ and F^y = C^y/F in B~~ — pp7r~ and B~ —> AAtt”, 
respectively. When t = (pe + PbY ^min = ^aa — ^pp /^- Furthermore, the 

relations from the numerator of the form factor are ~ (7^/2 and Cy^ ~ —Gy/lO. Once 
we combine the relations above and integrate them through the phase space, we obtain that 
Br{B~ —AAtt”) ~ O{10~^)Br{B~ — pp7r~). 

Similarly, we can also study the neutral decay modes of B^ —>■ AAK^{7r^), which have 
not been measured yet. The decay branching ratios are found to be 

^ AAK^) = (2.5 ± 0.3) x 10"®, 

Br{B^ AAtt®) = (0.4 ± 0.4) x 10"^ (23) 

As seen from Eq. (121, the decay branching ratio of B^ —> AAJF® is almost the same as that 
of B~ AAK~ , whereas B^ —»■ AAtt® is still suppressed as B~ —>■ AAtt”. We note that the 
errors of Br{B~ —> AAK~) and Br{B^ —> AAK^) in Eqs. (|22|l and (El are small since the 
main contributions are not from which receive almost all uncertainties from the data. 

To describe the possible non-factorizable effects, we also £t the data with N^. = 2 and oo. 
As expected, we End that the branching ratios of the K modes are slightly changed, while 
the central values of B~ —> AAvr” and .B® —> AAtt® shift to 2.3 and 1.3 (0.7 and 0.1) for 
Nf. = 2 (oo), respectively. We may conclude that the two tt modes remain small even with 
including all possible non-factorizable effects. However, as pointed in Ref. j^, 02 can only 
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be determined by the experimental data in the two-body B decays, since the experimental 
data of Br{B^ —>■ and Br{B^ are much larger than the theoretical values, 

which means the failure of the factorization approximation. On the other hand, in the three- 
body baryonic B decays, due to the complicated topology of Feynman diagrams, it is not as 
easy as those of the two-body decays to influence 02 by the annihilations j3l| as well as the 
hnal state interactions Therefore, the value of 02 may not change much. Nevertheless, 
we leave the surprise to the experimentalists if the factorization does not work well in these 
two modes as those in the two-body decays. 

In sum, we have studied three-body charmless baryonic B decays of i? —> AAvr and 
B —> AAA' based on the QCD counting rules in the standard model. We have shown that 
Br{B~’^ —>■ AAK~’^) S> Br{B~’^ —> AAtt^’^) unlike the cases of i? —*■ ppn (K). Explicitly, 
we have found that Br{B~ AAK~) = (2.8±0.2) x 10“®, Br{B^ AAK^) = (2.5±0.3) x 
10“®, Br{B- AAtt-) = (1.7±0.7) x 10“^ and Br{B° AAtt^) = (0.4 ±0.4) x 10“^ It is 
interesting to note that the decay of B^ —> AAK^ should be seen at the ongoing B factories 
soon. 
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